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Abstract 

The Left-Right symmetric extension of the Standard Model with Higgs isospin triplets 
can provide neutrino masses via a TeV scale seesaw mechanism. The doubly charged 
Higgs bosons and induce lepton flavour violating decays — > III at tree-level 
via a coupling which is related to the Maki-Nakagawa-Sakata matrix (I^vIns)- We study 
the magnitude and correlation of — > /// and ji ^ e^y with specific assumptions for the 
origin of the large mixing in Vmns while respecting the stringent bound for /i eee. It is 
also shown that an angular asymmetry for III is sensitive to the relative strength of 

the Hf^^ and mediated contributions and provides a means of distinguishing models 
with doubly charged Higgs bosons. 



PACS index: 13.35.-r, 14.60.Pq, 14.80.Cp 

Keywords : Higgs boson, Neutrino mass and mixing, Lepton flavour violation 



* akeroyd@inail.ncku . edu . tw 
'''mayumi. aoki@kek.jp 
yasuhiro . okada@kek.jp 



1 Introduction 



In recent years there has been increasing evidence that neutrinos oscillate and possess a small 
mass below the eV scale [1]. This revelation necessitates physics beyond the Standard Model 
(SM), which could manifest itself at the CERN Large Hadron Collider (LHC) and/or in low 
energy experiments which search for lepton flavour violation (LFV) [2]. Consequently, models 
of neutrino mass generation which can be probed at present and forthcoming experiments are 
of great phenomenological interest. 

Massive neutrinos may be accommodated by adding a SU{2)l singlet ("sterile") right- 
handed neutrino z/^ to the SM Lagrangian together with the corresponding Dirac mass term. 
In order to obtain masses of the eV scale, the Yukawa coupling of the neutrinos to the SM 
Higgs boson would need to be at least 6 orders of magnitude smaller than the electron Yukawa 
coupling. Moreover, there would be no observable phenomenological consequences aside from 
neutrino oscillations. More appealing frameworks for neutrino mass generation can be found if 
neutrinos are of the Majorana type. The celebrated seesaw mechanism [3] ascribes the smallness 
of the neutrino mass to the large scale of the unobserved heavy right-handed neutrinos {Nji). 
Dirac mass terms of the order of the top quark mass would require Nji ~ 10^^ GeV, a scale 
which is far beyond the reach of any envisioned collider. Reducing the scale of Nji to the 
order of a few TeV would require Dirac mass terms of the order of MeV, which constitutes 
a mild fine-tuning with respect to magnitude of the charged lepton masses. However, such a 
choice would permit the mechanism to be probed at future high-energy colliders. This "low 
energy seesaw mechanism" may be implemented in Left-Right (LR) symmetric models [1] with 
Higgs triplet representations [5], in which the mass matrix for is given by the product of a 
new Yukawa coupling hij and a triplet vacuum expectation value (vev) vr. The mass scale of 
the other new particles in the model (e.g. the new gauge bosons Wr, Zr and doubly charged 
Higgs bosons H^^) is also determined by vr, resulting in a rich phenomenology at the LHC if 
VR^TeV [6], [7].' 

The Yukawa coupling hij mediates many low energy LFV processes. In this paper we 
consider the impact of on the branching ratio (BR) of the LFV decays r hljlk and 
/X — >■ 67 in the context of the LR symmetric model [8j,[9j. Experimental prospects for fi ^ 
are bright with the imminent commencement of the MEG experiment which will probe BR~ 
10~^^ — > 10~^^, two to three orders of magnitude beyond the current upper limit [10]. At the 
e~^e~ B factories limits of the order BR(r — hljlk) < 10~^ with ~ 90 fb^^ |ll],[l2] have been 
obtained utilizing direct e+e~ r'^r^ production. Simulations of the detection prospects 
at a proposed high luminosity e+e~ B factory with £=5— > 50 ab~^ anticipate sensitivity to 
BR~ 10~^ 10~^ [13j. Additional searches can be performed at the LHC where r leptons are 
copiously produced from the decays of W, Z, B, D, with anticipated sensitivities to BR~ 10~® 

In the LR symmetric model mediate r Uljlk at tree-level due to an effective 4- Fermi 
charged lepton interaction proportional to h*ihjk/M^±±. Hence such a model can comfortably 
accommodate BRs of order 10~^ — >■ 10~^ which will be probed at current and forthcoming 
experiments. For the loop induced decays fi cy and t I'-f the dominant contribution 
in the LR symmetric model originates from diagrams involving H^^. In the LR model one 
has BR(r III) ^ BR(r Z7), which contrasts with the general expectation BR(r — >■ 



I'f) ^BR(r — s> ///) for models in which the tree-level r kljlk interaction is absent (for 
scenarios where BR(r — I'j) ~BR(r III) is possible see [15]). Due to the larger backgrounds 
for the search for r — ^ I'y, the experimental sensitivity to BR(r I'-f) at the e^e~ B factories is 
expected to be inferior to that for r kljlk pTSJ- Consequently the hierarchy BR(r III) ^ 
BR(r — >■ /7) in the LR model affords more promising detection prospects. 

The presence of the above tree-level 4-Fermi interaction would also mediate the decay /i — >■ 
eee for which there is a strict bound (< 10~^^ fT6]) at least three orders of magnitude stronger 
than the anticipated experimental sensitivity to r ^ III. Hence obtaining BR(r —>■ III) > 10~® 
together with compliance of the above bound on BR(/i eee) restricts the structure of hij. In 
the LR model we consider a specific ansatz for hij motivated by the observed pattern of the 
neutrino mixing angles, and perform a numerical study of the magnitude of BR(r — kljlk)- 

Observation of BR(r — > ///) > 10~^ would be a spectacular signal of LFV and could be 
readily accommodated by a tree-level 4-Fermi coupling such as /i*j/ijfc/M^±± in the LR model. 
However, other models which contain a or with an analogous h^j leptonic Yukawa 
coupling can cause a similar enhancement of BR(r III). If a signal were established for 
r — > /// the angular distribution of the leptons can act as powerful discriminator of the models 
[17]. Such studies could be carried out at a high luminosity e+e~ B factory [13j . 

Our work is organized as follows. In section 2 the manifest LR symmetric model is briefly 
reviewed. In section 3 a numerical analysis of BR(r — > ///) and BR(/i 67) is presented. In 
section 4 we discuss how to discriminate the LR model from other models which contain a 
by means of angular asymmetries in the LFV decays. Conclusions are contained in section 5. 

2 Left-Right Symmetric Model 

The Left-Right (LR) symmetric model is an extension of the Standard Model (SM) based on 
the gauge group SU{2)fi (g) SU{2)l ® U{1)b-l- The LR symmetric model has many virtues, 
e.g. i) the restoration of parity as an original symmetry of the Lagrangian which is broken 
spontaneously by a Higgs vev, and ii) the replacement of the arbitrary SM hypercharge Y by 
the theoretically more attractive B — L. Although the Higgs sector is arbitrary, a theoretically 
and phenomenologically appealing way to break the SU{2)ji gauge symmetry is by invoking 
Higgs isospin triplet representations. Such a choice conveniently allows the implementation 
of a low energy seesaw mechanism for neutrino masses. The vev of the neutral member of 
the right-handed triplet [vji) can be chosen to give a TeV scale Majorana mass term for the 
right-handed neutrinos, while the bidoublet Higgs fields provide the small Dirac mass, leading 
to light masses for the observed neutrinos. The above LR symmetric model predicts several 
new particles, among which the new gauge bosons Wr, and doubly charged scalars -ff^^, 
[ig have impressive discovery potential at hadron colliders if f r = (9 (1-10) TeV due to 
their large cross-sections and/or low background signatures. 

Experiments which search for LFV decays of /i and r provide a complementary way of 
probing the LR symmetric model. A comprehensive study of /x — 67, /i eee and \i ^ e 
conversion in the present model was performed in [9]. However, the recent termination of 
the MECO {fi e conversion) experiment together with no immediate improvement for the 
SINDRUM collaboration limit BR(/i eee) < 10"^^ [T6] leaves /i — > 67 as the only means of 
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testing the LR model in LFV processes involving fx in the near future [10] . 

An alternative probe of the LR model which was not developed in [9] are the LFV decays 
r hljlk- Although the experimental sensitivity is inferior to that for the above processes 
involving /x, the decays r hljlk have the virtue of probing many combinations of the triplet 
Higgs-lepton-lepton Yukawa hij couplings. We introduce various strutures for the arbitrary 
hij motivated by the currently preferred bi-large mixing form of the Maki-Nakag aw a- S akat a 
[T9] (MNS) matrix. Should a signal for r — >■ III and/or /i ^ 67 be observed, the angular 
distribution of the final state leptons can provide a means of distinguishing models with if^^. 
We now briefly introduce the LR model and present the relevant formulae for the numerical 
discussion in Section 3. For a detailed introduction we refer the reader to [20] . 

The quarks and leptons are assigned to multiplets with quantum numbers (T^, Tr, B — L): 

Q.L = f j 1 : (1/2 : : 1/3) , Q,^ = ^| ^ : (0 : 1/2 : 1/3) , 





:(l/2:0:-l), Un = \ ) : (0 : 1/2 : -1) . (1) 

R 

Here % = 1,2,3 denote generation number. The spontaneous symmetry breaking to U{l)em 
occurs through the Higgs mechanism. The Higgs sector consists of a bidoublet Higgs field, 
and two triplet Higgs fields, Al and Ar: 

The vevs for these fields follows: 



The gauge groups SU{2)ji and U{1)b-l are spontaneously broken at the scale vr. Phenomeno- 
logical considerations require k = y^^f + ~ Hil^fL (EW scale). The vev vl does not 

play a role in the breaking of the gauge symmetries and is constrained to be small {v^ < S 
GeV) in order to comply with the measurement of p = Mz cosOw / My/ ~ 1. The LR model 
predicts six neutral Higgs bosons, two singly charged Higgs bosons, and two doubly charged 
Higgs bosons. The Lagrangian is required to be the invariant under the discrete left-right sym- 
metry: Ql ^ Qr , Ll ^ Lr , Al Ar , $ ^ $"1". This ensures equal gauge couphngs 
{Ql = 9r = 9) for SU{2)l and SU{2)r. 

The leptonic Yukawa interactions follows: 

- ^J'"^- = LLivD^ + yumR + lyuiLlCT^ALLL + LICt^ArLr) + h.c. . (4) 

Here $ = r2$*r2; t/d, yn are Dirac type Yukawa coupling; um is a 3 x 3 Majorana type Yukawa 
coupling matrix which will lead to Majorana neutrino masses (see below) and is the primary 
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motivation for introducing the Higgs triplet representations A^, and Ar. Invariance under the 
left-right discrete symmetry gives Vd = Vd, Vd = Vd and t/M = Vm- Redefinitions of the fields 
Ll and Lr enable yM to be taken as real, positive and diagonal, while maintaining yo = y\), 
Vd = Vd- Hereafter yM is taken in this diagonal basis. The 3x3 mass matrix for charged 
leptons is: 

Ml = -^{yDK2 + Vdi^i) , (5) 

which is diagonalized by the unitary matrices, ¥[ and V^, as: 

V^j^MiVli = diag{me, m^, m^) . (6) 
The Lagrangian for the neutrino masses is: 

- Cn^ass = ^{nLM^riR + fiRMlriL) , (7) 

where ul = {i^l, ^r)^ and hr = (z/£, ur)"^ with the definition of z/|j = C{urY . The 6x6 mass 
matrix for the neutrinos can be written in the block form: 



M.=[ r/ )■ (8) 

V vn\y Mr ' ^ ' 



Each entry is given by: 



"^D = (yDfi^i + yDi^2) , (9) 

Mr = V2yMVR, (10) 
Ml = V2yMVL • (11) 

The neutrino mass matrix is diagonalized by a 6 x 6 unitary matrix V as V'^ MjjV = M^*"^ = 
diag{mi,m2, m^, Mi, M2, M3), where rrij and Mj are the masses for neutrino mass eigenstates: 

The small neutrino masses rrii are generated by the Type II seesaw mechanism. Obtaining 
eV scale neutrino masses with yM = 0(0.1 — 1) requires Ml (and consequently vl) to be eV 
scale. However, the minimization of the Higgs potential leads to a relationship among the vevs, 
Vl ~ 'Jk^/vr, where 7 is a function (introduced in [7]) of scalar quartic couplings Pi and pi. 
For natural values of Pi and pi one has 7 ~ 1 and thus vl would be 0(1 — 10) GeV for vr ~ 
TeV. Reducing fL to the eV scale to order to comply with the observed neutrino mass scale 
would require severe fine-tuning 7 < 10^^. In LR model phenomenology it is standard to set 
Pi = (and hence 7 = 0) which ensures vl = 0. Henceforth we will take vl = for which 
the masses of the light neutrinos arise from Type I seesaw mechanism and are approximately 
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rrii ~ m'jj/Mji. In order to realize the low energy (~ (9(1 — 10) TeV) scale for the right- 
handed Majorana neutrinos, the Dirac mass term ttid should be O (MeV), which for K2 ~ 
corresponds to yo ~ 10^^ (i.e. comparable in magnitude to the electron Yukawa coupling). 

There are two physical singly charged Higgs bosons, and H^, which are linear combi- 
nations of the singly charged scalar fields residing in ^,Ai and A/j. The leptonic couplings tjo 
of (which is essentially composed of 0f and 0f ) are of order mi/m\y and can be neglected 
compared to leptonic Yukawa couplings for the triplet field ~ 6^ which are unrelated to 
fermion masses and may be sizeable. The interaction of iff with leptons is as follows (where 
Nl = V^riL, Nr = V'^nji, N = Nl + Nji = N'^ and I = Ir + II are the neutrino and charged 
lepton fields respectively in the mass eigenstate basis, and Pl,r = (1 =F 75)/2): 



H+ N (h Pj) l + H{ I (h^ Pr) N 



(13) 



The LFV interactions of leptons with doubly charged Higgs bosons (where = 5^^, 
Sr"^ for Vl = 0) are given by: 



Ctt±± 



Ht% I' ihL,R Pl,r) I + H-^;r I {h^R Pr,l) I' 



The LFV coupling matrices in Eq. f[T3l) and Eq. ffT^ are respectively given by: 

1 / vr^ ' 



h 

hL 
hR 



V2VR 

1 

V^Vr 

1 

V^Vr 



L 



MrVI, 



v'^MrvI, 



V'^MrVr . 



(14) 



(15) 
(16) 
(17) 



Note that /i is a 6 x 3 matrix and and Hr are 3x3 matrices. 
The mass matrix for the charged vector bosons is: 



-2kiK2 + 2vr 



:i8) 



This is diagonalized via the mixing angle ^ = —tan ^ {2k,ik,2/vr) /2 with the eigenvalues 



= 9' + vlT ^4 + ^^14) /4: 

Wl = cos ^Wi + sin ^ W2, 



W, 



R 



sin ^Wi + cos ^ W2 . 



(19) 



The strong experimental constraint on the mixing angle (^ < 10 ^) [31] enforces one of k,i,k,2 
to be small if vr = O (TeV). Neglecting such small mixing between Wi and W2, the LFV 



^ For an alternative approach which maintains the Type II seesaw mechanism and obtains ~eV by means 
of horizontal symmetries see |21) . 
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interactions with the gauge bosons follows: 



cc 



i-wt, + i 



7" Pr (kI) 



N-W. 



2 n 



+N 



'r Pl {Kl)] I ■ Wt^ + Ik Pl {kI)] N ■ W{ 



(20) 



where Kl and Kr are the 6x3 LFV coupling matrices which are respectively written as: 



K, 



^MNS 



^R 



(21) 



The upper 3x3 block in Kl can be identified as the hermitian conjugate of the MNS matrix 
Vmns on neglecting C^(^) contributions. 

2.1 Manifest LR Symmetric Model 

In the LR model the mixing matrices for the left and right fermions are in general not equal 
e.g. for the lepton sector Vl ^ V^. The special case of Vl = is referred to as the "Manifest 
LR symmetric model" and arises in either of the following scenarios: i) both ni and k,2 are 
real, or ii) one of ki and K2 is identically zero. In our numerical analysis we will set k,2 = 0, 
which has the virtue of eliminating Wl — Wr mixing and in some cases (for specific forms of 
the Higgs potential) is required to suppress FCNCs and preserve unitarity in the LR model [6]. 
In the Manifest LR symmetric model one has the additional constraint niD = rn^L) which must 
be respected when evaluating the magnitude of the LFV processes. 

A further important consequence of the Manifest LR symmetric model is the relationship 
fiL = liR = h which can be derived from Eq. (fT6i) and Eq. lfTTI) . Using Kl and Kr in Eq.(l2Tl). 
the LFV couplings for the interactions of leptons with singly and doubly charged Higgs bosons 
are as follows: 

h = K*h , h= -^—Kl Mf""^ Kr . (22) 



At leading order in itld/Mr, one may express h by: 

n=heavy 



Xr, 



\V2vr) 



(23) 
(24) 



2.2 Effective Lagrangian and branching ratios for the LFV processes 
2.2.1 4-lepton interactions 

The effective Lagrangian for 4-lepton interactions is as follows: 



'^ = 7;ih*)rm{h)jk 



M' 



\- (jml'-PLlk) {hl^^PLlj) + 



H 



Hi 



(jn^rPRh) (kl^PRlj 



(25) 
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The branching ratio for r Uljlk is given by: 

85 



BR{t -> kljlki 



h^,? ( + 1 BR{t ^ /iH 



(26) 



Here 5*=! (2) for j = k {j 7^ fc). In the LR model one may express hij in terms of K^K^ via 
Eq. (l23|) . However, the above form Eq. fl26|l can be apphed to other models with Hf^^ and i/^^ 
for which no identity Eq. (|23|1 exists. 



2.2.2 /i^e7 

The effective Lagrangian for — > 67 is as follows: 

AG 



^ {m^An-pa^"'PLeF^^ + m^ALW" PrcF^^ + h.c.} 



Al receives contributions from W2 — Ni and and is given by [9]: 

r2 _ T\/r2 



where 



A, 



, , ,Ki] (Kf 

n=heavy 



E 



2~ '-'Sl^^rij 



W2 



3 M^, 



2;(l + 2x) ^ 



'1 — X + logx) + 1 



8(1-X2)2 4(1 -X)2 [(1-X)2 

Api^ receives contributions from Hf;^ and Hf and is given explicitly by [9j: 



A 



R 



1 



n=heavy 

The branching ratio for /i 67 is given by: 

5/?(/i ^ 67) = 3847rV(|Ai|2 + |ArJ2) 



24 Ml± 



(27) 



(28) 



(29) 



(30) 



(31) 



3 Numerical analysis for BR(r III) and BR(/i 67) 

The LEV decays r — >■ ZZZ are the analogy of fi —>■ eee and provide sensitive probes of the h^j 
couplings in the LR model. Mere observation of such a decay would constitute a spectacular 
signal of physics beyond the SM. There are six distinct decays for r+ III (likewise for r~): 
r+ n'^ , r+ e~^e'^e~, r+ — /i+/i+e^, r+ ^ yU+/i~e+, r+ — > e~*"e~'"/i~, r+ e^e^fi^. 
Searches for all six decays have been performed by BABAR (91 fb~^) [12j and BELLE (87 fb~^) 
[TT] . Upper limits of the order BR(r III) < 2 x 10~^ were derived. Although these limits 
are several orders of magnitude weaker than the bound BR(/i — > eee) < 10"^^, they have the 

^ In our numerical analysis we do not include a suppression factor of ^ 15% arising from electromagnetic 
corrections 
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virtue of constraining many combinations of the hij couplings in the context of the LR modeL 
Moreover, greater sensitivity to BR(r III) is expected from forthcoming experiments. A 
proposed Super B Factory anticipates sensitivity to BR(r III) ~ 10^^ and 10~^ for 5 ab^^ 
and 50 ab~^ respectively [13] • At the LHC, r can be copiously produced from several sources 
(from B/D decay and direct production via pp ^ W ^ tu, pp ^ Z ^ r+r") and sensitivity 
to BR(r — i> ///) > 10~* is claimed [14J. Such low BRs can be reached due to the very small SM 
background. In contrast, the background to r — > Z7 is non-negligible and might prevent the B 
factories from probing below BR~ 10^®. In addition, one expects BR(r — > III) ^ BR(r — > Ij) 
in the LR symmetric model and thus the former decay is the more effective probe. We note 
that other rare LFV r decays involving quark final states will not arise in the LR model since 
hij mediates processes involving leptons only. Hence we shall only focus on r — ///. In our 
numerical analysis only the stringent constraint from /i — eee is imposed. Other constraints 
on hij (e.g. the anomalous magnetic moment of yU {g ~ 2), Bhabha scattering and other LFV 
processes - see [22]) are considerably weaker and are neglected. 

The magnitude of hij cannot be predicted from the neutrino oscillation data alone since it 
is related to the physics at SU(2)j:j breaking scale. However, h^j also crucially depends on the 
mixing matrix in the charged lepton sector ¥[: 

^ VI^MrVI. (32) 



V2vR 

We have neglected the C'(^) contribution and used the convention that Mr is diagonal and 
positive. Since Vl also enters the MNS matrix: 

V^MNS = VjJV^ , (33) 

we will introduce 4 distinct structures for ¥[ motivated by the bi-large mixing form of Vmns 
and perform a quantitative analysis of the magnitude of hij (and consequently BR(r — ^ ///)) in 
the LR model. 

In order to establish our formalism we first explicitly present the standard parametrisation 
of the MNS matrix: 



MNS 



/ ci3 sise 
1 
-sne'^ ci3 



/ 1 

C23 S23 

V -S23 C23 

U{92■,)U{9n)Ui9^2) • (34) 




Here Sij (cjj) represents sin6'jj {cos 9ij) and the unitary matrices U{92z), U{9i^), and U{9i2) are 
responsible for mixing between 2-3, 1-3, and 1-2 elements respectively. 

The angles 9i2 and 6*23 are measured with relatively good accuracy in the solar and at- 
mospheric neutrino oscillation experiments respectively. The solar and KamLAND reactor 
neutrino oscillation experiments [23], [21] provide the following constraints on the mixing angle 
9i2 and the mass-squared difference of Am^g — T^\—ni\'- sin^ 9i2 ~ 0.31 , ^rn^^ ~ 8x 10^^ eV^ . 



The mixing angle ^23 and the mass-squared difference Am^g measured in the atmospheric neu- 

-^'131 



trino oscillation are as follows [25], [26]: sin^ 2^23 ~ 1-0 , lAmfgl ~ 2.6 x 10"^ eV^ . An upper 
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bound on the remaining angle 613 has been obtained from the CHOOZ and Palo Verde reactor 
neutrino oscillation experiments [27], [28]: sin 6*13 ^0.2 . 

The ignorance of the sign of Am^g and the absolute neutrino mass scale leads to the fol- 
lowing three neutrino mass patterns which are consistent with current oscillation data: Normal 
hierarchy (NH) mi < m2 <^ m^; Inverted hierarchy (IH) -C mi < 7712] Quasi degeneracy 
(DG) mi ~ 1712 ~ iTT-s- Data from WMAP [29j provides the following constraint on the sum of 
the light neutrino masses: J2i=i,2,3 rrii < 2 eV. However, LFV processes in the LR model do not 
depend sensitively on the neutrino mass pattern. 

In order to perform our numerical analysis of the magnitude of BR(r III) and BR(/i — 67) 
we introduce the following four specific cases: 



CASE 






I 




ll 


II 


-7U{923)U{9is) 


iu{ei2) 


III 


~iu{e23) 


iu{9i3)u{ei2) 


IV 


-il 





Here / represents a unit matrix. In CASE I (IV) both large mixings in Vmns originate from the 
charged lepton (neutrino) mixing matrix. Each case has distinct ways of satisfying the stringent 
bound BR(yU — > eee) < 10^^^. In our numerical analysis we will assume multi-TeV scale masses 
for Hf^^ and -ff^^ which renders direct detection improbable at the LHC. For M^±±,M^±± < 1 

TeV the LHC has excellent discovery prospects in the channels H^^ if if [18] (especially 
for lij = e,yu). Observation of (which would provide a measurement of Mh±±) together 
with a signal for r — > kljlk would permit a measurement of the coupling combination 

We wish to study the magnitude of BR(r III) and BR(/i 67) in the parameter space 
with a phenomenologically acceptable neutrino mass matrix. In the LR model the light neutrino 
masses arise from the seesaw mechanism and are approximately as follows: 

= —m£)M]^^mJ) . (35) 

At the leading order, rriu is diagonalized by V^: 

m^ ~ Vi;mf''3Vl^ , (36) 

where mf'°'^ = diag{mi, 7712, m^). The Dirac mass matrix mo depends on an arbitrary Yukawa 
coupling yo- As advocated in |30], it is beneficial to parametrize a general seesaw type matrix 
such that the arbitrary m^ is replaced by potential observables i.e. the heavy and light neutrino 
masses. We will apply the formalism of [30] in our numerical analysis, with the additional 
constraint that the manifest LR model requires the Dirac mass matrix for both the neutrinos 
and charged leptons to be a hermitian matrix: 

mo = m^o ■ (37) 

We introduce a complex orthogonal matrix R which satisfies R^R = 1 and parametrize the 
neutrino Dirac mass matrix m^ as follows: 

mn = -iV^^mf'^'R'^R . (38) 
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The LFV processes are evaluated in the parameter space where there exists an R matrix which 
satisfies the condition Eg. (1371) . This condition guarantees a phenomenologically acceptable 
neutrino mass matrix and perturbative Yukawa coupling yo. In our calculation, the Majorana 
phases in the MNS matrix are neglected while the CP conserving cases for the Dirac phase 
(5 = or tt) are taken into account for simplicity. We will comment on the case of the CP 
violating Dirac phase in section 3.5. Neglecting CP violation, the condition Eg. fl37l) reguires 
that is purely imaginary, while i? is a real matrix. We stress that the LFV processes do not 
depend on the actual structure of R. However, proving the existence of an R matrix for each 
of the four cases (I, II, III, IV) ensures the validity of our numerical analysis. 

3.1 Numerical results: CASE I 

The bi-large mixing originates from the charged lepton sector. We parametrize the R matrix 
as follows: 

R = uie^,)uie^^uie^,) . (39) 

The explicit form in CASE I is given by: 



nR nR nR 

^12 ~ ^23 ~ ^13 



0, 



(40) 



which leads to R = I. In this case mo = y mi^"'^ Mr, and thus mo = rn^D is automatically 
satisfied. 




(a) 



10' 



BR(x"^|a laV") ^ 10 



BR(x^^e e^e^) 



-14 -13 -12 

10 10 10 10 



10 



10"'^^ 



10"' r 



10"" r 



(b): 

BR(x^^|l e^e^) -_ 



BR(x e |i ) 



10"'"^ 10"^^ 10"'^ 10 
BR(|a"'^e"'Y) 



10" 



10" 



10" 



10" 



10"'^ 



(c) 



BR(x^^e l^V^) 



BR(x^^|i e^|X^) 



10"'^ lO"^""* 10"'^ 10" 
BR(|^"'^e"'y) 



BR(|a y) 

Figure 1: Branching ratios for (a) r+ e^e^e^ and r+ —>■ fi'^^fi^, (b) r+ —>■ e~e^fi^ and 
r"*" /i"e"'"e+, and (c) r+ e~/i+/i+ and r+ yU^e+/i+ against BR(/i+ — > e'^'y) for S13 = 
in CASE I. The experimental bound on BR(/i eee) (< 10~^^) is imposed. 



Using Eg. (l32|) the elements of hij in CASE I are as follows: 

-V2vnhef, = Ci2Si2C23Ci3(M2 - Mi) ± Ci3Si3S23(M3 - - cj^M,) , (41) 
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Figure 2: Branching ratios for (a) r~*" ^ e e^e^ and r+ ^ /i fi~^fi^, (b) r+ ^ e e~^fi'^ 
and r+ — > fi~e~^e~^, and (c) r+ e^fi^fi^ and r+ — > /i^e^/i"*" against BR(/i+ e'^l) for 
sin = 0.2 with 5 = in CASE I. 



(M2 - Ml) ± C23C13S13 (M3 - - C22M1) 



-V^t^flV = C23S23(C?3M3 - 4M2 - S?2^i) ± Ci2Si2Si3((s^3 - 4)M2 + M^] 

+ C23S23Si3(Si2M2 - C^aMi) , 
-^VRKe = C?3(4Mi + sl^M^) + S^gMg , 

-V2VRh^^ = cj^isl^Mi + cl^M2) + sl^cl^Ms T 2CuSl2C23S23Sl3{M2 - Ml) 



(42) 

(43) 
(44) 

(45) 



Here the upper (lower) sign is for Dirac phase S = (vr). It is clear that the off-diagonal 
elements in the h couplings vanish when all the heavy neutrinos are degenerate in mass, i.e. 
Ml = M2 = M3. The strict bound on BR(/i — >• eee) requires |/ie^| -C 1 which leads to the 
following conditions: 



(i) Mi~M2, Si3(M3-M2) ^0 , 

(ii) S13M3 ~ Ci2Si2(M2 - Ml) + 513(4^2 + ^Mi) 



TT . 



(46) 
(47) 



In (i) both terms contributing to /ig^ are zero, while in (ii) there a cancellation which ensures 

We show the branching ratios for all six r"*" —>■ III decays against BR(/i"'" —>■ e'^'j) for 
Si3 = in Fig{Tl For simplicity we assume degeneracy for the masses of the heavy particles: 
= M^±± = M^±± = Mjj± = 3 TeV. In our numerical analysis we vary the heavy 
neutrino masses randomly in the range 1 TeV < Mj < 5 TeV, with distributions that are flat 
on a logarithmic scale. This range is consistent with the vacuum stability condition for vr given 
in [33] . 

In Fig{T](a) ((b), (c)) the light and the dark points respectively denote the branching ratios 
for r+ — >• e~e'^e~^ and — > ( r+ — >• e~e~^fi'^ and t'^ — >• fi~e~^e~^, r+ fi~e~^n~^ and 
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Figure 3: Branching ratios for (a) r~*" ^ e e~^e~^ and ^ fi fi~^fi^, (b) r+ ^ e e~^fi'^ 
and r+ — > fi~e~^e~^, and (c) r+ e^fi^fi^ and r+ — > /i^e^/i"*" against BR(/i+ ^ e"''7) for 
sin = 0.2 with 5 = vr in CASE I. 



r"*" — > e~^'^fi'^). We impose the experimental constraint BR(/i — eee) < 10~^^ which prevents 
a large mass difference between Mi and M2 when S13 = 0, as shown in condition (i) of Eq. p6|) . 
Among the six r"*" — >■ III decay modes the branching ratios of r+ — > fi~ fi'^fi'^ and r+ /i~e^e+ 
can reach the anticipated sensitivity (> 5 x 10^^) of a future B factory. Such branching ratios 
are realized when the mass splitting M3 — M2 assumes larger values because h^^r increases with 
it. On the other hand, the r — ^ e transition is suppressed because \he^\ = |/ier| for S13 = 0, and 
\hen\ is necessarily small in order to comply with the severe constraint from fi eee. This can 
be seen for the light points in Fig{T] (a) where BR(r — >■ eee) is proportional to BR(/x — eee). 
Moreover, there is a strong correlation BR(r eee) ~ 10 x BR(/i e'-f). BR(/i e'y) can be 
large as 10~^^, which is within the sensitivity of MEG experiment. 

Figs J2] and [3] show the branching ratios for r III for S13 = 0.2 with the Dirac phase 6 = 
and TT respectively, imposing the constraint BR(/i eee) < 10"^^. For the other parameters 
we take the same values as in FiglH When 6 = (Figl2]) all LEV processes are predicted to 
be small because of the small mass differences between the heavy neutrinos as shown in the 
condition (i). However, there is still the possibility of observing yu ^ e7 at MEG experiment. 
On the other hand, taking 6 = tt (FiglS]) results in observable branching ratios for r — > In 
this scenario one has |/ie^i| < l^erl- Therefore the equality of \he^\ and l/ierl in Fig {T] is broken, 
which enables enhancement of BR(r"'" e~e~^e^) and BR(r+ e^/i^/i"'') with simultaneous 
suppression of the /x — e transition. Moreover, BR(/i — e'j) can be large, resulting in multiple 
signals of LFV processes. 

3.2 Numerical results: CASE II 

In this case the large mixing for the atmospheric angle originates from the charged lepton 
sector, while the large solar angle originates from the neutrino sector. The R matrix (given in 
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Eq. (!39|) ) which satisfies the condition mn = rrif) is: 



tan = rjj ^ ^^ ^=tangi2 , (4^ 



el = ef3 = 0. (49) 



The exphcit form for /ijj is obtained from Eg. (1411) - fl45l) by taking Si2 = 0. The condition for 
suppressing \hf,^\ is: 

Si3(M3-4Mi)^0, (50) 

which requires S13 ~ or Mi ~ M2 ~ M3. In the latter case none of the r — > III decays are 
measurable, as in CASE I with S13 = 0.2 and 5 = 0. On the other hand, /ig^ and her are 
zero when S13 = 0, which results in vanishing branching ratios for r+ e~ l~^ l~^ , fi~ and 
yU — > 67. FigJH shows the branching ratios for r+ —>■ fi~fi^fi~^ and r'^ —>■ fi~e^e'^ against the 
heaviest neutrino mass M3. Clearly the branching ratios increase with M3 and for M3 > 3 TeV 
observable rates are attained. 




, : • ^R(T e e )- 

1 2 3 4 5 

M (TeV) 



Figure 4: Branching ratios for r+ ^ /i^yU^ and ^ jj e^e^ against M3 for sin 6*13 = in 
CASE II. 

When Si3 7^ the magnitudes of each entry in hij are the same for S = and vr. We show 
in Figl5]the branching ratios for r III decay against BR(/i — 67) for S13 = 0.2. In order to 
satisfy the condition in Eg. (1501) there cannot be large splittings among the masses of the heavy 
neutrinos. In this scenario only BR(yU 67) reaches the future experimental sensitivity. 



3.3 Numerical Results: CASE III 

In CASE III the constraint from fi eee is satisfied automatically since /le^ = (obtained by 
setting S12 = Si3 = in Eq. fHTl) . We obtain R as follows by treating 6^13 as a perturbation: 

tan9« = ^^±^tan«.,, (51) 
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Figure 5: Branching ratios for (a) r+ ^ e e+e+ and r+ — /i'^/U^, (b) r+ e e^^^ and 
r+ — > /i^e+e"*", and (c) r+ jj,^ ji'^ and r+ — > ji^e^n^ against BR(/i — > 67) for S13 = 0.2 



with (5 = 0, TT in CASE II. 



sin 6'§ = f si3 



/ 



(52) 



+ V'MiM2m3 , (53) 
5f = \JlvU,Ci2Si2 I (^\J M2mim2 + ^iW^mi^ cfg^ - M2mim2 + ^^^7712"] Sj 



'12 



+ |y^M2M3(m2S?2 - "^icD + \jMiMzmim2[ 



■^12 '-I2J r '-12'='12 



+v^1 \^MiM2m2 + MsJ^] S12C 



'12 



MiM2mi + M3 C12S 



'12 r ' 



sin^f^. 



13 " 



^M3(Y/mrCi2sf2 — V^Si2cf2)s^3 + {v/M3^+ V^( ■\/^Ci2cf2 + ^/rn^Si2sf2) \ Sl3 



VMiTUs + VM3(ym7Ci2Cf2 + A/^'5i2Sf2 



(54) 



.(55) 



CASE III predicts the same results as CASE II for the branching ratios of the LEV processes 
with si3 = (EiglH). 



3.4 Numerical results: CASE IV 

In this case the bi-large mixing originates from the neutrino sector. We checked the existence 
of R numerically. It is clear that none of the LEV processes are observed since hij is a diagonal 
matrix. 
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Table 1: LFV processes within the sensitivity of forthcoming or planned experiments for CASES 



I, II, III and IV with Mw^ = M^±± = Mj^±± = Mj^± = 3 TeV. BR(/i ^ 67) > IQ-^^ and 
BR(r ^ > IQ-^ is denoted by V". 



CASE 


sin 013 
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fi'^ — !■ e+7 















V 
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V 
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0.2, 
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0.2, 


TT 
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V V 
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II 
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TT 
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III 


0.2, 
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0.2, 


TT 
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IV 


0.2, 















0.2, 


TT 











3.5 Summary 

The numerical results of CASES I, II, III and IV for Mw^ = iVf^±± = M^±± = M^± = 3 TeV 
are qualitatively summarized in Table 1. It is clear that the number of observable rates for the 
LFV decays r —>■ III and /i — > 67 depends sensitively on the origin of the bi-large mixing in 
with up to five signals being possible in the optimum scenario (CASE I with sin ^^13 = 0.2, 
6 = 7c). Hence future searches for r III and fi —>■ e'j provide insight into the structure of hij 
in the LR model. 

Both BR(r —>■ III) and BR(/i — > eee) are inversely proportional to the fourth power of Mh±± 
(when Mh±± = M^±± = M^±±) as shown in Eg. (126!) . Reducing Mw2i Mjj±±, and M^±± in 

L R L Ft 

all the figures would cause enhancement of BR(r —>■ III) and BR(/i — > eee) while maintaining 
the correlation. However, BR(yU — > eee) vanishes in CASE II with S13 = and CASE HI, and 
hence the strong bound BR(/i —>■ eee) < 10~^^ is automatically satisfied for any Mu±±. In these 
latter cases the current experimental limits for BR(r — > III) can be reached and thus an upper 
bound on the heaviest neutrino mass M3 can be derived {e.g. FiglHwith = Mh±± = 1.5 
TeV gives BR(r ^ /i/i/i) > 10"^ for M3 > 1.8 TeV). 

We note here that the presented numerical results are for the scenario of all phases in 
^MNS taken to be zero. We now briefiy discuss the effect of including a non-zero Dirac phase 
{5 7^ 0) in our analysis El- In CASE I and CASE II the presence of 5 7^ would not affect 
Eq. fl38|) since U{6i^) in Eq. flM|) (which contains 5) only appears in Vl and not in V^. Thus 
the solution for the R matrix in the CP conserving case (5 = or tt) also holds for the 

The presence of Majorana phases in ^mns would effectively change the relative phase of the heavy neutrino 
masses Mi in Eqs.(41) to (45). In such a scenario small \hee\ can arise from a cancellation among the terms 
in Eq.(44), which provides an additional way to suppress BR(/x — s- eee) in CASE I. We performed an explicit 
numerical analysis and found a pattern of lepton flavour violation different from the case of |/ie/i| ~ 0. In 
particular, BR(t+ e~e^/x+) and BR(r+ ji^e^ji^) can be enhanced to observable rates. We thank the 
referee for bringing this scenario to our attention. 
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case of 5 7^ 0. Since Vl (which determines the LFV rates) has some dependence on 6 we 
would expect some changes in our numerical results. For example, we have calculated the 
branching ratios for the LFV processes for 6 = 7i/2 and 3tt/2 with S13 = 0.2 in CASE I 
and II. For all cases BR(yU 67) < 2 x 10^^^. Maximum values of 10~^ were found for 
BR(r+ ^ /i-/iV^), BR(r+ ^ /i-e+e+), and BR(r+ ^ e'i2+fi+), with smaller (< 10"^°) 
values for BR(r+ e~e+e+), BR(r+ e~e+/i+) and BR(r+ — > /i"e+/i+). When 6^0 there 
is another possibility to suppress BR(yU — > eee) by cancellation in h^e, because the sign of the 
term sf^M^ in Eg. (1441) flips for 6 = 7t/2 and 3tt/2. However, the cancellation is not so significant 
because of the small factor 5^3. In CASE III and IV the presence of 5 7^ increases the number 
of free parameters in and no general solution for the R matrix can be found. However we 
expect that Eg. ( 1371) will be satisfied in specific regions in the parameter space of Mi, and for 
CASE HI some observable LFV rates might still be possible. 



4 P odd Asymmetry for r kljh cind n ^ 

From the preceding sections it is evident that both BR(r III) and BR(;U 67) can be 
enhanced to experimental observability in the LR symmetric model. In particular, BR(r 
III) > 10~^ would be a signal suggestive of models which can mediate the decay at tree-level 
e.g. the LR model via virtual exchange of 0. In order to compare the LR model with 
other models we introduce the Higgs Triplet Model and Zee-Babu Model in section 4.1 and 4.2. 
These models provide a different mechanism for neutrino mass generation, and can enhance 
both r hljh and /i — > 67 to experimental observability by virtual exchange of or 
H^^. If a signal were established for any of the six decays r — > III, further information on the 
underlying model can be obtained by studying the angular distribution of the leptons. In section 
4.3 we show how the P odd asymmetry for both decays may act as a powerful discriminator of 
the three models under consideration. 



4.1 Higgs Triplet Model 

In the Higgs Triplet Model (HTM) [35] a single / = 1, F = 2 complex SU{2)l triplet is added 
to the SM. No right-handed neutrino is introduced, and the light neutrinos receive a Majorana 
mass proportional to the left-handed triplet vev {vl) leading to the following neutrino mass 
matrix: 



Mu = V2vl 



h^,e hf,f, h^r I • (56) 

y h-fe hrfi hr-j 

In the HTM hij is directly related to the neutrino masses and mixing angles as follows: 

hij = -^^V^^^diagimi, 1712, m3)V^^^ ■ (57) 

Formally, this expression for hij is eguivalent to that in CASE I with the replacements (mi, 7712, m^) 
(Ml, M2, M3) and vl — vr. In Eg. (1571) vl is a free parameter and is necessarily non-zero (un- 
like in the LR symmetric model) in order to generate neutrino masses. Its magnitude may 



For models with large extra dimensions see 
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lie anywhere in the range eV < < 8 GeV, where the lower limit arises from the require- 
ment of a perturbative hij satisfying Eq. (!57|) and the upper limit is derived from maintaining 
p ~ 1. In the HTM there is no H^"^ and so r — Uljlk is mediated solely by Hf^^. Obtaining 
BR(r -H> liljlk) > 10~* with mu±± = 1 TeV requires \h*^hjjs\ > 10"^. Due to the ignorance of 
vl the magnitude of hij cannot be predicted and so the HTM can only accommodate observable 
BRs for T liljlk- However, unlike the LR symmetric model, the HTM provides predictions 
for the ratios of r — >■ Uljlk [36] with are distinct for each of the various neutrino mass patterns 
NH, IH and DG. The necessary suppression of BR(/i — > eee) relative to BR(r III) requires 
hen sufficiently small. As in CASE I Eq.( l47|) . this is arranged by invoking a cancellation 

between two terms contributing to /le/^, one depending on 6*13 and the other depending on both 
9i2 and r = Am^g/^^ia [36],[37j. Observation of r — >■ /// would restrict ^13 to a narrow interval 
which can be predicted in terms of 9i2 and r. For the case of NH neutrinos ^13 ~ \/r while for 
IH and DG neutrinos 6*13 ~ r. Since current oscillation data suggests r ^ 0.04, the value of 6*13 
need to ensure small hp., is of the order 0.1 for NH and 0.01 for IH and DG. 



4.2 Two loop radiative singlet Higgs model (Zee-Babu model) 

Neutrino mass may be absent at the tree-level but is generated radiatively via higher order 
diagrams involving L = 2 scalars. In the Zee-Babu model (ZBM) SU{2)l singlet charged 
scalars H^^ and are added to the SM Lagrangian [5B],[SS] with the following Yukawa 
couplings: 

= k {l^lCL%) eatHt + {iIrCI.r) + h.c. (58) 

No right-handed neutrino is introduced. A Majorana mass for the light neutrinos arises at the 
two loop level in which the lepton number violating trilinear coupling jj.H^H^H'^'^ plays a 
crucial role. The explicit form for A^^, is as follows: 



/ 2 \ 



\ 



ee'iu, 



e/i 



(59) 



where e = feTlfiir-, ^ = fe/i/ fur, ^ij = hijUiirrij {mi.rrij are charged fermion masses), hij 
KjC^Kj) for ^ = i 7^ i) and C is given by: 



C 



(60) 



Here J is a dimensionless quantity of 0{1) originating from the loop integration. Clearly the 
expression for Ai^ differs from that in the HTM and involves 9 arbitrary couplings. Since the 
model predicts one massless neutrino (at the two- loop level), quasi-degenerate neutrinos are not 
permitted (unlike the HTM) and only NH and IH mass patterns can be accommodated. The 
/ couplings (contained in e and e') are directly related to the elements of A4^. In the scenario 
of NH, e ~ e' ~ tan6'i2/v^ and sin ^13 is close to zero. Since e, e' < 1 one may neglect those 
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terms in Ai^ which are proportional to the electron mass (i.e. Uee, ^^e/i, i^er)- This simplification 
leads to the following prediction [39]: h^^ : h^r '■ hrr ~ 1 : rn^/nir : {m^/nir)^. In the case of 
IH, large values are required for e, e'(> 5), and thus neglecting tUee, i^e/x, i^er in A^^ may not be 
entirely justified. However, if such terms are neglected then the above prediction for the 
ratio of h^^ : h^^ : h^r also approximately holds for the case of IH. In the ZBM there is no 
and so r — i> hljlk is mediated solely by H^^. 

Another significant difference with the HTM is that eV scale neutrino masses requires /, 
h^^ ~ 10~^, and thus LFV decays cannot be suppressed arbitrarily if the 2-loop diagram is 
solely responsible for the generation of the neutrino mass matrix. Such relatively large couplings 
are necessary since a rough upper bound on ( (which is a function of model parameters) can 
be derived. In contrast, in the HTM is arbitrary and eV scale neutrino masses can be 
accommodated even with hij ~ 10~^°. The requirement that /, /i^^ ~ 10~^ suggests that 
BR(/i — i> 67) and BR(r — > /x/x/x) could be within range of upcoming experiments [lO] . 

Since h^e, h^^ and her may be treated as free parameters (essentially unrelated to the 
neutrino mass matrix) the necessary suppression of /i — ^ eee can be obtained by merely choosing 

and/or /ige very small. Observable rates for BR(r — eee) can be arranged by appropriate 
choice of hee and her- 



4.3 Sensitivity of P odd asymmetry to Hf^ and 

Angular distributions of LFV decays can act as a powerful discriminator of models of new 
physics. The predictions for /i+ — Cj^-^ and jX^ e^7 depend on the chirality structure of 
LFV interactions and so in general would be model dependent. Ref. |41j defined various P 
odd and T odd asymmetries for /i — e7 and fi eee and performed a numerical analysis in 
the context of supersymmetric SU{5) and SO{10). Analogous asymmetries were defined for 

— > /^7 and — ^ /// in [T7j. In this section we apply the general formulae introduced in 
Refs.|17].|41j to the three models of interest which all contain H^^. 

For the decay yU^ —>■ e~e^e^ with polarized fi~^, one defines 6^- as the angle between the 
polarization vector of /i+ and the direction of the e~, the latter taken to be the z direction. 
The P odd asymmetry Ap is defined as an asymmetry in the 6'e- distribution. In contrast, for 
T produced in the process e"'"e~ — *• t~^t~ the helicity of the r in the LFV decay r —>■ III is not 
known initially. Consequently, the experimental set up is sensitive to both III or —>■ III. 

However, by exploiting the spin correlation of the pair produced r (i.e. e+e~ '^l'^r^'^r'^l) 
information on the helicity of the LFV decaying r can be obtained by studying the angular 
and kinematical distributions of the non-LFV decay of the other r in the r — > event. For 
illustration we shall always take the non-LFV decay mode as r — vrz/, although such an analysis 
can also be performed for other main decay modes such as r — pz/, Oiz/, IvV. 

In the notation of [TTj the effective 4-Fermi interaction for r+ fi^fi^fi^ mediated by 
is as follows: 

-4Gf 

= {g;iT-i^PRii){jIl^.PRl^) + g^iT-i^PLli)ijIl^.PLli)} . (61) 

The differential cross-section for the events r+ p^/i+p+ (LFV) and — tt^z/ (non-LFV) 
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Table 2: Expressions for g^, in the three models 



HTM(gr) ZBM(gg") LR(ggi) 



4Gf „ 













V2 ^4 











KP 



is: 



= a[e e t t )BR[t ^ -k u) I ^28"^ j — 2 — dcosO 

S — 2?Tl'^ 

X ^{Fcos^lcos^^ 

s + 2m:^ 



(62) 



where 



X = (Ifi-aP + \94\'^)ai{xi,X2); Y = (l^igp - \g4\'^)ai{xi,X2) . (63) 

Here ai{xi,X2) is a function of the energy variables Xi = 2i?i/m^ and X2 = 2E2/mr where 
Ei[E2) is the energy of ^'^ with larger (smaller) energy in the rest frame of r"*": 

ai{xi, X2) = 8(2 - xi- X2){xi + X2 - 1) . (64) 

The angles 9 and (f) specify the decay plane of r+ ji^ n^jji^ relative to the production plane 
of e^e~ T^T~ in the r"*" rest frame. The angle 6'^ is the angle between the direction of 
momentum of r~ and vr" in the r~ rest frame. For a detailed discussion we refer the reader to 
[T7] . It is clear that Y determines the angular dependence of Eq. and thus F is a measure 
of the P odd asymmetry for r /i/i/i. For our quantitative study of its magnitude we define: 

A{t ^/i/i) = / ^ . (65) 

Clearly A{t — >■ /i/i/i) = (±1) corresponds to zero (maximal) asymmetry. 

The expressions for g^ and g^ in the three models under consideration are given in Table 2. 
In the manifest LR symmetric model hij cancels out in the expression for A{t —>■ III) leaving a 
simple dependence on Mjj±± and M^±± which applies to all six r III decays: 

L R 

l/MU± - l/MU± 



In FiglH] (a) A{t —>■ III) in plotted in the plane {Mjj±±, M^±±). Clearly the case of degeneracy 
{Mjj±± = Mjj±±) gives A{t III) = 0, while M^±± > Mjj±± {Mjj±± < M^±±) results in 

L Ft L R L R 

positive (negative) A{t — > ///). For the HTM and ZBM the asymmetry is maximal, being 
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Figure 6: (a) A{t III) and (b) BR(r — > /i/x/i) and in the plane {Mj^±±, 



±± 



— 1 and +1 respectively. Consequently, A{t — * ///) in the LR symmetric model may differ 
significantly from the corresponding value in models with only a (HTM) or i?^^ (ZBM). 
Thus A{t — > III) has the potential to discriminate between the various models with a if^^. 
In addition, in the context of the LR model a measurement of A{t — > ///) provides important 
information on the ratio Mjj±±/M^±±, which could assist the direct searches for H"^^ and i/^^ 
at high energy colliders. For comparison we show in Figl6] (b) the dependence of BR(r /i/i/i) 
on M^±± and M^±± for = 0.05. 

Of the order of 50 r ^ /// events would be needed to distinguish A{t III) = +1 from — 1. 
A high luminosity upgrade of the existing B factories anticipates up 10^*^ r+r^ pairs and thus 
BR(r — s> > 10~^ would allow measurements of A{t — > ///). From Figl6] (a) it is clear that 
a LR model with M^±± <^ Mjj±± {Mjj±± -C Mj^±±) will give an almost maximal A{t III) 
and consequently would be difficult to distinguish from the HTM (ZBM). However, if a signal 
were also observed for /x+ the analogous P odd asymmetry, A{fi ej), can serve as 

an additional discriminator: 

In CASE I, Figj3] there is a parameter space for BR(r nnn) ^ 10~* and BR(/i 67) ~ 
10~^^), which might provide sufficient events for both asymmetries to be measured. In contrast 
to r ^ /^yU/^, the loop induced decay — > 67 can be mediated by Hf^, Hf and in the 
LR symmetric model. One may write a simplified formula for Al and Ar (written explicitly 
in Eqs. ( |28l) and ( |30l) ). where a,b,c,d are functions of masses and the heavy neutrino mixing 



matrix Kj^: 

Al = a (M^±± ,KR)+b (M^± , Kn, M,) , 

An = c (M^±± ,KR) + d {Mh± , Kr) . (68) 



In the LR model usually a,c^ b,d and so the dominant contribution to A{fi 67) arises from 
H^^ and H^^. Hence in LR models one expects A{fx 67) ~ A{t III). This is shown in 
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Table 3: P odd asymmetries A{t —>■ III) and A{fi ej) in the three models 







HTM (H^^) 


ZBM (Hp) 






67) 


-1 


-1 < ^ < +1 


-1 < ^ < +1 


A{r- 


///) 


-1 


+ 1 


-1 < ^ < +1 



FigI7]where A{t nnn) is plotted against Ai^ji 67) for CASE I with sin 6'i3 = 0.2 and 5 = -n. 
The plotted points correspond to observable rates for both LFV decays (taken to be 10^^ < 
BR(r /i/i/i) < 10~^ and 10~^^ < BR(/i 67) < 10~^^), within the following parameter 
region: Mw^ = 3 TeV, 1 TeV < < 5 TeV, 2 TeV < M^±± = M^± ^ M^±± < 4 TeV. 
Each point also satisfies the constraint BR(/i eee) < 10^^^. Clearly the vast majority of the 
points are close to the line A{ij 67) = A{t III), showing that the diagrams involving 
and H^^ give the dominant contribution over most of the parameter space. The asymmetries 
differ sizeably only when Mj and Mw^ are considerably smaller than M^±± and M^±± . In the 
HTM one has: 

= 0, 

Ar = c{Mj,±±,h,,)+d {M^±,hij) . (69) 
and thus A{fi 67) is maximal. In the ZBM: 

Al = a {Mfj±±,hij) , 

Ar = diM^lj,,) . (70) 

A{fi 67) may take any value since the masses of and are unrelated and hij 7^ fij 
in general. The allowed ranges of A{t III) and A{fi 67) in the three models under 
consideration are summarized in Table 3. It is clear that if signals for both r III and fi ^ e'j 
are observed, the corresponding asymmetries may act as a powerful discriminator of the models. 



5 Conclusions 

The Left-Right symmetric extension of the Standard Model with TeV scale breaking of SU (2)^ 
via a right handed Higgs isospin triplet vacuum expectation value provides an attractive ex- 
planation for neutrino masses via the seesaw mechanism. The doubly charged scalars 
and with mass of order TeV mediate the LFV decays r — /// at tree-level via a Yukawa 
coupling hij which is related to the Maki-Nakagawa-Sakata matrix (Vmns)- We introduced 
four ansatz for the origin of the bi-large mixing in Vmns which satisfy the stringent bound 
BR(/i — i> eee) < 10"^^ in distinct ways. A numerical study of the magnitude and correlation of 
BR(r^ — > ///) and BR(yU 67) was performed. It was shown that the number of observable 
rates for such LFV decays depends sensitively on the origin of the bi-large mixing in I^^ns, 
with multiple LFV signals being possible in specific cases. 

If a signal for r —>■ III were observed we showed how the definition of an angular asymmetry 
provides information on the relative strength of the contributions from and H^^. Such an 
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Figure 7: Correlation between A{f^ 67) and A{t /^/^A*) in the LR model. 

asymmetry may also be used to distinguish the LR symmetric model from other models which 
contain either or and thus predict maximal asymmetries. 
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